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ABSTRACT 

We present a photometric study of the globular cluster systems of the Fornax cluster galaxies 
NGC 1374, NGC 1379, and NGC 1387. The data consists of images from the wide-field MO- 
SAIC Imager of the CTIO 4-m telescope, obtained with Washington C and Kron-Cousins R 
filters. The images cover a field of 36 x 36 arcmin, corresponding to 200 x 200 kpc at the 
Fornax distance. Two of the galaxies, NGC 1374 and NGC 1379, are low-luminosity ellipti- 
cals while NGC 1387 is a low-luminosity lenticular. Their cluster systems are still embedded 
in the cluster system of NGC 1399. Therefore the use of a large field is crucial and some dif- 
ferences to previous work can be explained by this. The colour distributions of all globular 
cluster systems are bimodal. NGC 1387 presents a particularly distinct separation between 
red and blue clusters and a overproportionally large population of red clusters. The radial dis- 
tribution is different for blue and red clusters, red clusters being more concentrated towards 
the respective galaxies. The different colour and radial distributions point to the existence of 
two globular cluster subpopulations in these galaxies. Specific frequencies are in the range 
Sm = 1-4 — 2.4, smaller than the typical values for elliptical galaxies. These galaxies might 
have suffered tidal stripping of blue globular clusters by NGC 1399. 

Key words: galaxies: individual: NGC 1374, NGC 1379, NGC 1387 - galaxies: clusters: 
general - galaxies: elliptical and lenticular, cD - galaxies: star clusters - galaxies: photometry 
- galaxies: haloes 



1 INTRODUCTION 

There is an extensive literature on globular cluster systems (GCSs) 
of lar ge galaxies loc ated in the ce ntre of clus ters (e.g. the reviews 
from lAshman & Zep f 1998 or (Harrij|200lll but low-luminosity 
galaxies with their own GCSs are also present in such environ- 
ments. The centres of galaxy clusters are interesting regions for 
studying possible interactions between the central massive galaxy 
and its less massive galactic neighbours. These interactions may 
also influence the properties of their cluster systems. 

The nearby Fornax cluster offers the best opportunity for 
studying such a scenario due to the huge GCS that s urrounds its 
central giant galaxy, NGC 1399 toirsch et al]l2003a '. and refer- 
ences therein), and the presence of several low-luminosity galax- 
ies in its neighbourhood. The GCS of NGC 1404, one of thes e less 
massive galaxies, has been studied b y iRichtler et alj il992h and , 
more recently, bv lForbes et alj il99d) and Grillmai r et alj Jl999h . 
The latter authors identified two GC subpopulations: metal-poor 
GCs (blues) and metal-rich ones (reds), the metal-rich clusters 
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being more centrally concentrated than the metal-poor ones. The 
low specific frequency Sn (as defined by Harris & van den Ber^ 
1981) of this system has been understood as a deficiency of GCs, 
probably stripped in an interaction process with NGC 1399 (Forbes, 
Brodi e & Grillmair 1 997). Numerical simulations were performed 
bvBe kki etalj<2003h to analyse this scenario of tidal stripping and 
accretion. 

One may speculate if GCs lost by low-luminosity galaxies are 
really captured by the central galaxy or just remain unbound within 
the galaxy cluster , like the prop osed intracluster GCs (White 198^ 
iWest et alil995l:lBassinol;tall20()3lh . 

Our aim with the present paper is to study more early-type 
galaxies in the Fornax cluster: NGC 1387 (type SO), NGC 1379 
(EO) and NGC 1374 (El). These galaxies are located at angular dis- 
tances from NGC 1399 of 19, 29 and 41 arcmin, respectively. 

The GCS of NGC 1387 has firs t been studied with CCD tech- 
niques bv 'Greb el & Richtled jl990^ who obtained a surprisingly 
red mean colour {B — V ) = 1.0. Later, Kohle et al. (1996) and 
iKissler-Patig et all |l997a) performed V and I photometry of its 
GCS as well as of NGC 1374 and NGC 1379. The former paper 
focused on the GC luminosity functions, the latter on the colour, 
spatial distributions, and total numbers of GCs. They concluded 
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that the galaxies have between 300 and 500 GCs (5'jv « 4 ± 1) 
and that their {V — I) colours are redder than those of the Milky 
Way GCs. However, they found no evidence for multiple popula- 
tions, but simply pointed to a wide range in metallicities. None of 
the GCSs seem to be elongated and the density profiles follow the 
respective galaxy light profiles. The background was determined 
locally within the frames of 7.4 x 7.4 arcmin used in these papers. 

The GCS of NGC 1379 has been the target of a HST-study 
iElson et'al]ll998h combining HST data (B band) with ground- 
based observations in the B and I bands. They estimated the to- 
tal population in NGC 1379 to be 440 GCs, based on a Gaussian 
luminosity function with mean < _B > = 24.95 and as = 1-55 
mag. Again, no evidence for bimodality was found in the (B — I) 
colour distribution. A specific frequency 5'jv ~ 3.5 was obtained 
up to a radius of 70 arcsec . The background corrections were mea- 
sured in a field located 1?4 from the Fornax cluster centre (see also 
[Grillmair et al. 1999). 

The apparent lack o f multiple GC subpopulations in these 
three galaxies was taken bv lKissler-Patig. Forbes & Minnitilil998h 
as an argument to set several constraints on the process of forma- 
tion of elliptical galaxies. The impact posed by the existence of GCs 
subpopulations on the studies of origin and evolution of galaxies is 
one of the reasons why we think that their GCSs deserve to be re- 
examined with a better metallicity res olution. 

On the basis of data pr o vided by iKissler-Patig et all ^9973), 
iForbes. Brodie & GrillmaiJ il997h suggest that the early-type 
galaxies in Fornax, located close to NGC 1399, might have been 
tidally stripped off their GCs. They note that their Sn values in- 
crease with projected distance from the cluster centre. Moreover, 
their GC surface density slopes are similar to that of the underly- 
ing starlight, perhaps as a result of losing the outermost clusters. 
Both effects have also been studied through dynamical models by 
iBekkietalJilOol . 

On the other hand.^ Capuzzo-Dolcetta & Donnar umma' ^20011) 
use the Elson et al. ( 1998) data for NGC 1379 and the Forbes etal] 
Il998 ) data for NGC 1399 and NGC 1404 (all from HST) to show 
that the surface brightness profiles of the underlying galaxies are 
steeper than the GC profiles. They propose that GCs from the in- 
ner regions might have been lost as a consequence of dynamical 
friction and tidal interaction with the galactic nuclei. 

We re-investigate the GCSs in NGC 1387, NGC 1379 and 
NGC 1374, on the basis of new MOSAIC images in the Washington 
photometric system. This system provides a metallicity resolution 
sufficiently good for a firm detection of GCs subpopulations. The 
large field improves the determination of the background (partly 
due to the NGC 1399 GCs) that may otherwise alter the radial den- 
sity profiles and luminosity functions of these low number GCSs. 

The analysis of the GCSs of our three target galaxies is or- 
ganized as follows: Section 2 describes the observations and the 
adopted criteria for the globular cluster candidates' selection. In 
Section 3 we analyse the colour and spatial distributions, the lu- 
minosity functions of the GCSs and also estimate the total popula- 
tions. A discussion of the results is provided in Section 4. Finally, a 
summary of the results and their implications are provided in Sec- 
tion 5. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Observations and point sources selection 

The observations were performed with the MOSAIC camera (8 
CCDs mosaic imager) mounted at the prime focus of the 4-m 
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Figure I. C + R combined image of the MOSAIC field. Labels indicate 
the three target galaxies. North is up and East to the left. 



Blanco telescope at the Cerro Tololo Inter-American Observatory 
(CTIO), during 2001 November 17-19. One pixel of the MOSAIC 
wide-field camera subtends 0.27 arcsec on the sky, which results in 
a field of 36 x 36 arcmin (200 x 200 kpc at the Fornax distance). 
All three target galaxies were covered by one MOSAIC field. For 
more information on the MOSAIC camera we refer to its homepage 
I http://www.noao.edu/kpno/mosaic/mosaic.htmI I. 

Kron-Cousins R and Washington C filters were used. We 
selected the R filter instead of the original Washington Ti as 
iGeisIeJ il996h has shown that the Kron-Cousins R filter is more 
efficient than Ti, and that R and Ti magnitudes are very simi- 
lar, with just a very small colour term and zero-point difference 
(R-Ti^ -0.02). 

In order to fill in the gaps between the 8 individual MOSAIC 
chips, the data were dithered taking three images in R with ex- 
posure times of 600 s each, and six images in C with exposures 
of 1 200 s each. Fig. 1 shows a C + R combined image of the 
MOSAIC field where the three galaxies have been identified. Their 
basic data are given in Table I. 

The MOSAIC data were reduced using the mscred package 
within IRAF^. The software corrects for the variable pixel scale 
across the CCD which might cause a 4 per cent variability in point 
sources brightness (from the centre to the comers). The final R 
image resulted with remaining sensitivity variations up to 0.8 per 
cent, and the final C image up to 1.5 per cent. The seeing in these 
final images was 1 arcsec on the R frame, and 1.2 arcsec on the C 
frame. 

The halo light of each galaxy was subtracted by means of a 
ring median filter with an inner boundary of 1.3 arcsec and an outer 
boun dary of 10.7 arcsec. I n previous work based on MOSAIC im- 
ages jPirsch et alj2003ilbh . we have already checked that this pro- 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which is operated by AURA, Inc. under contract to the National Science 
Foundation. 
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Table 1. General data for the target galaxies in the 
Fornax cluster. Total blue apparent magnitudes (Bt) 
and heliocentric radial velocities (RV) were taken 
from the RC3 catalog, and redden ing values from 
ISchlegel. Finkbeiner & DavisHl998ft . 



Galaxy 


RA(J2000) 


DEC(J2()00) 


Type 




RV[kmsec-M 


E{B ~ V) 


NGC 1374 


3'"35™17'' 


-35°13'35" 


El 


12.00 


1351 ± 15 


0.014 


NGC 1379 


3'"36™03'' 


-35°26'28" 


EO 


11.80 


1376 ± 19 


0.012 


NGC 1387 


3''36"57'' 


-35°30'24" 


SO 


11.68 


1328 ± 24 


0.013 



cess does not affect the photometry. Anyway, we also have com- 
pared the magnitudes obtained with and without filtering for several 
point sources on different positions on the image. The differences 
were negligible as compared to the photometric errors (see below). 

The photometry was done with DAOPHOT within IRAF. For 
the first search we used DAOFIND on a combined C + -R image 
and 20 856 unresolved and extended objects were detected in the 
whole field. In the final C and R images, a second order variable 
PSF was derived using about 100 evenly distributed stars per frame, 
which was adjusted to the sources through the ALLSTAR task. The 
estimated errors in the aperture corrections between the PSF radius 
and a 15pix radius were ~ 0.02 and 0.01 mag for the C and R 
images, respectively. 

Several tests were carried out before deciding how to select the 
point sources. We compared the results obtaine d using the stellar- 
ity in dex estimated by the SExtractor software iBertin & AmoutsI 
Il996h which is a reliable tool for separating point so urces and 
extended objects (see, for instance, iNonino et alJll99^ . adopting 
an stellarity in dex bet ween 0.4 and 1 for point sources (follow- 
ine lDirsch et al ."2003aX with the results obtained from DAOPHOT 
using different options for the parameters x ^nd sharpness of the 
ALLSTAR task. We finally adopted as point sources those objects 
selected with the ALLSTAR parameters \ < 2 and sharpness be- 
tween -0.5 and 0.5 (see, for instance, Reid et al. 1996). 

In this way, out of the total number of objects originally de- 
tected in the field, 6 520 point sources were selected, populating a 
magnitude range of 18.5 < Ti < 24.5. 



2.2 Photometric calibration and globular clusters selection 

The photometric calibration of our observations already has been 
performed byjjirsch et al. 1 2003s) so we conveniently can use their 
relations. Fiel ds of standard stars were selected from the list of 
iGeisleil il996h . Four or five fields, with roughly 10 standard stars 
located on each of them, were obtained in each night. Care was 
taken to cover a large range in airmasses (1.0 to 1.9 approximately). 

In order to calculate the standard magnitudes from the in- 
strumental ones, we have used the equations already obtained by 
lDirschetalJ<2003ah for this run, which we reproduce here: 



0.3 
0.2 
0.1 
0.0 



21 22 
TI 



Figure 2. Photometric eiTors in T\ from DaoPhot plotted against the T\ 
magnitude for all point sources in the MOSAIC field. 
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Ti = Rinst + (0.72 ± 0.01) - (0.08 ± 0.01)Xr 
+ (0.021 ± 0.004)(C-Ti) 
C = Cinst + (0.06 ±0.02)- (0.30 ±0.01)Xc 
+ (0.074 ± 0.004) (C-Ti) 

These equations also include the zero-point difference to 
transform R into Ti magnitudes (as explained in Section 2.1) so, 
in the rest of this paper, we will use Ti magnitudes, not 7?. The 
photometric errors for the point sources are shown in Fig. 2. 

The colour-magnitude diagram for all point sources in the field 
is plotted in Fig. 3. Globular cluster candidates can be clearly dis- 
tinguished in the range 0.8 ^ (C — Ti) ^ 2.3 mag. Objects bluer 
than (C — Ti) = 0.8 and fainter than Ti = 23 are mostly compact 



C-T1 

Figure 3. Colour-magnitude diagram for all point sources in the MOSAIC 
field. Globular cluster candidates show up within the box defined by the 
colour range 0.8 ^ (C — Ti) ^ 2.3 and magnitude range 20 ^ Ti ^ 
23.7. 



background galaxies while objects redder than (C — Ti) ~ 2.8 are 
mainly foreground stars. 

We select as GC candidates those point sources with colours 
0.8 < (C - Ti) ^ 2.3 and magnitudes 20 ^ Ti 23.7 where 
the faint Ti limit is adopted according to the completeness re- 
sults described in Section 2.3. All GC candidates identified on the 
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Figure 4. Projected distribution of GC candidates on the MOSAIC field: 
the colour {C — T\) = 1.45 is adopted as limit between blues (open cir- 
cles) and reds (crosses). The GCSs around the three target galaxies can be 
seen as well as the background distribution decreasing from East to West 
that corresponds to the NGC 1399 GCS. Solid lines indicate the regions se- 
lected as background for each galaxy (labelled 'a', 'b' and 'c') and large 
circles the regions excluded from them (see text). The MOSAIC scale is 
0.27 arcsec pixel" ^ . 



the region (labelled 'c' in Fig. 4) located at y < 3000 pix and 
X < 6000 pix (density 1.42 ± 0.06 objects arcmin"^) where we 
excluded the region with y < 3000 and x > 6000 pix due to the 
presence of NGC 1373, a small galaxy at 5 arcmin to the NW with 
respect to NGC 1374 (see Fig. 1), and whose probably own cluster 
candidates show up at a: = 8000 and y = 1000 in Fig. 4. The back- 
ground defined for NGC 1379 is the region (labelled 'b') located at 
3000 < y < 5000 pix excluding a circle of radius 800 pix (density 
1.53 ± 0.07 objects arcmin"^); for NGC 1387 the region (labelled 
'a') located at a; < 7000 and y > 5000 pix excluding a circle of ra- 
dius 1000 pix (density 2.79 ± 0.09 objects arcmin"^) was selected 
as background; the region with x > 7000 and y > 5000 pix was 
excluded due to the presence of another galaxy, NGC 1380B. 

All these regions are shown in Fig. 4. It can be seen that the 
background density increases from West to East as we get closer to 
NGC 1399. This effect is also present in the V background counts 
obtained for the same galaxies by Kohle et al. 1 1996, see their 
fig. 1), although the influence of NGC 1399 GCs is not mentioned 
as a background component. 

The reddening towards our three galaxies taken from 
ISchlegeletai] <1998^. is depicted in Table 1. The colour excess in 
the Washington photometric system E{C' — Ti) is estimated by the 
relation E{C - Ti) = 1.97E{B - V) (Harris & Canterna 197^ 
and th e absorption in R with Ar/Av ~ 0.75 I Rieke & Lebof skM 
1 198^. For the reddening of the three background regions defined 
in the previous paragraph, we adopt the values corresponding to 
the galaxies located in each of them, respectively. In the rest of the 
paper, all magnitudes and colours will be reddening free. 



Table 2. Catalog of globular cluster candidates identified on the MOSAIC 
field. The full version of this table is available in the online article. 



RA='(J2000) 


Dec.''(J2000) 


Ti 


(C-Ti) 


3''34'"38!78 


-35°19'18'.'6 


23.02 ± 0.05 


2.22 ±0.09 


3 34 39.55 


-35 17 27.4 


21.78 ± 0.03 


2.11 ±0.04 


3 34 40.64 


-35 23 01.7 


21.52 ± 0.02 


1.37 ±0.02 



^ The estimated precision is 0.3 arcsec. 



MOSAIC field are shown in Fig. 4, where the GCSs around the 
three galaxies can clearly be seen; their positions, magnitudes and 
colours are listed in Table 2. We also notice in Fig. 4 a smooth 
background distribution of GC candidates with a projected density 
decreasing from East to West, representing the GCS of NGC 1399, 
the central galaxy of the Fornax cluster. It is located at a projected 
angular distance of 19 arcmin from NGC 1387, towards the East. 

The selection of GC candidates should be corrected for the 
contamination by the background. In the study of the NGC 1399 
GCS we have used for this purpose a back ground fie ld located as far 
as 3? 5 from the parent galaxy IDirsch et £1112003 j) . In the present 
case, where there is a rather strong contamination from NGC 1399 
itself (see Fig. 4) we must estimate the background contribution in 
regions closer to the target galaxies, i.e. within the same MOSAIC 
image. We should take into account that blue clusters from the GCS 
of NGC 1399 extend up to 40 - 45 arcmin from the Fornax cluster 
centre (Bassino et al. 2005, in preparation) and that NGC 1374, the 
target galaxy located farther away from it, has an angular separation 
of 41 arcmin. Thus, we have chosen as background for NGC 1374 



2.3 Completeness 

In order to estimate the completeness of the GC candidates for the 
R and C final images, we added 1 000 artificial stars by means 
of the ADDSTAR task within DAOPHOT, which were generated 
with random x, y positions scattered through the whole MOSAIC 
field. They were distributed within approximate the same the colour 
range used to select the GC candidates. This process was repeated 
10 times, with 10 different seed numbers. The whole reduction pro- 
cess then was applied to each frame with added stars in exactly the 
same way as to the science images. The differences between in- 
put and output magnitudes and colours for the 'added' stars did not 
show any systematic trend with the magnitude. 

As a result, the completeness for stars in the colour range cor- 
responding to the GC candidates is close to 90 per cent for a limit- 
ing magnitude Ti — 23 decreases rapidly to 50 per cent for a lim- 
iting magnitude Ti = 23.7; this latter magnitude has been adopted 
as the faint magnitude limit of the GC candidates for further study 
in this paper. No difference was found in the completeness function 
for blue and red GCs, taking (C — Ti) — 1.45 as the colour limit 
between them. 

The analysis of the regions close to the galaxies' centres shows 
that inside r = 100 pix. (27 arcsec) the detection of GC candidates 
is affected by the light of the galaxies and saturation so we will not 
consider those inner regions in our study. The areas corresponding 
to 100 pix < r < 200 pix (27 - 54 arcsec) are still affected by 
the light of the galaxies and we expect there a completeness factor 
around 80 ± 5 per cent. The GC radial distributions will be cor- 
rected for this last effect and, in the following, we will refer to it as 
the innermost annulus or bin 'geometric incompleteness'. 



Globular clusters in Fornax Cluster galaxies 5 



3 THE GLOBULAR CLUSTER SYSTEMS 

In what follows we will analyse the GCSs of the three galaxies. On 
the basis of the projected distributions depicted in Fig. 4 and the ra- 
dial density profiles (see below), the radial extensions of the GCSs 
will be taken as r = 3.1 arcmin (TOOpix) for NGC 1387, r = 2.7 
arcmin (600 pix) for NGC 1379, and r = 2.2 arcmin (500 pix) for 
NGC 1374. We adopt an inner limit of at r = 27 arcsec (100 pix) 
for all GCSs (see Section 2.3). The edge-on disc galaxy NGC 1375 
is seen in projection close to NGC 1374 so its GCs may be con- 
taminating the NGC 1374 GCS; this is not causing a problem be- 
cause the NGC 1375 GCS does not seem to be very populous, as 
proved by inspection of the images. NGC 1373 seems to have its 
own cluster system, as mentioned above, so we restrict the selec- 
tion of GC candidates in NGC 1374 to a radial distance of 500 pix 
to exclude it. Table 3 gives the number of blue and red clus- 
ter candidates found in the regions defined for each galaxy. The 
colour range is 0.8 ^ (C — Ti) ^ 2.3 and the magnitude range 
20 ^ Ti 23.7, The selected colour limit between blue and red 
clusters is (C — Ti) — 1.45, as supported by the colour distribu- 
tions discussed below (Fig. 6). We remind here that the adopted 
faint magnitude cut-off is Ti = 23.7, the magnitude at which the 
completeness function is 50 per cent. 



3.1 Colour distributions 

Fig. 5 presents the colour-magnitude diagrams for point sources 
(selected as explained in Section 2.1) around the three galaxies 
within the radial limits stated above. GC candidates show up in 
them: the presence of two GC populations is clear in NGC 1387 
and NGC 1374. The separation of blue and red GC candidates is 
not so clear in NGC 1379. 

We can get a better picture from the colour distributions shown 
in Fig. 6. Displayed are the raw colour distributions, the back- 
grounds, and the background corrected distribution, all smoothed 
by a Gaussian kernel with dispersion 0.09 mag. The histograms of 
the background-corrected data are also plotted. Not only the GCSs 
of NGC 1387 and NGC 1374 appear clearly bimodal, but the one 
of NGC 1379 too. It also can be seen in Fig. 6 that the contam- 
ination by the background decreases with increasing angular dis- 
tance from NGC 1399 (angular distances of NGC 1387, NGC 1379 
and NGC 1374 from NGC 1399 are 19, 29 and 41 arcmin, respec- 
tively) and that this effect is mostly due to a blue population. This is 
in agreement with the less concentrated radial distribution of blue 
clusters found at large distances from the parent galaxy (Bassino 
et al. 2005, in preparation) with respect to the radial distribution 
of red GCs. In addition, we note that the colour distribution of the 
background selected for NGC 1399 by Dirsch et al. (20033, their 
Fig. 6), which is a field located 3? 5 northeast of the Fornax cen- 
tre, looks very similar to the background chosen for NGC 1374, 
the most distant galaxy from NGC 1399. The backgrounds for 
NGC 1379 and NGC 1387 present larger blue populations, demon- 
strating the large extent of the NGC 1399 system. 

Apparently, the better metallicity resolution provided by the 
Washington photometric system allows us to detect bimodali- 
ties that were not visible in other photometric systems. The 
GCSs of NGC 1387, NGC 1379 and NGC 1374 did not show ev- 
idence f or multiple popul ations in the {V — /) photometry by 
fKissler- Patig et al]il997a|), neither did the {B — I) photometry of 
NGC 1379 presented bv lsison et alj il998h . as already mentioned 
in the introduction. We are aware of the low metallici ty resolution 
of the (V — /) colour index (see, for instance, fig. 7 in lPirsch et alj 
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Figure 5. Colour-magnitude diagrams for point sources ai'ound the galaxies 
indicated in each figure. Vertical lines indicate mean colours corresponding 
to the blue and red GC candidates, respectively, taken from Table 3. 



l2003j) but one may wonder why the (B—I) data, known to provide 
a better metallicity resolution than {V — /) colours, do not show 
any evidence of bimodality for NGC 1379 either. In order to analyse 
the colour d istributions, the ba ckground corrections play a funda- 
mental role. lF,lsonetalj r998h used for their HST observations a 
background field lying 1?4 to the south of the centre of the Fornax 
cluster; it seems that this field probably lacks the strong contribu- 
tion added to the background by the NGC 1399 GCS. Also, in the 
same paper, there is no background field for the CTIO (ground- 
based) observations, so the authors use a corrected sample of the 
CFRS (Canada-France Redshift Survey, Lillv et al. 1995) which 
might probably present the same problem. 

In addition to the individual GCS colour distribu- 
tions of NGC 1374, NG C 1379, NGC 1387 and NGC 1427, 
iKissler-Patig et alJ h997iJ) presented,a {V — I) composite his- 
togram for all GCs with colour uncertainties less than 0.1 mag. 
Interestingly, this composite distribution seems to be bimodal 
according to a KMM test (the GCS of NGC 1427 has already 
been found to be clearly bimodal on the ba sis of the Washington 
photometry performed bv Forte et alJi200ll) ). 

The GCSs of NGC 1387 and NGC 1374 are more populous 
than the one of NGC 1379 within the quoted ranges of magnitude, 
colour, and radial distance (it should be noticed that vertical scales 
are different in Fig. 6). The number of identified red cluster candi- 
dates in NGC 1387 outnumber the blue ones by more than a factor 
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Table 3. Number of identified GC candidates (0.8 ^ (C — Ti) ^ 2.3 and 20 Ti ^ 23.7), and results of the Gaussian fits to their background-corrected 
colour distributions. 



Galaxy Nj^f^^ (C-Ti)biuc 'tmuc N^^^ N^°- ^ (C-Ti),,d <y.oA 

NGC 1374 75 ± 9 61 ± 9 1.27 ±0.02 0.11 ±0.02 81 ± 9 73 ± 9 1.70 ±0.03 0.19 ±0.03 

NGC 1379 64 ± 8 40 ± 8 1.28 ± 0.02 0.15 ±0.02 59 ± 8 49 ± 7 1.65 ± 0.01 0.10 ±0.01 

NGC 1387 94 ±10 32 ± 10 1.25 ± 0.01 0.08 ± 0.01 123 ± 11 100 ±11 1.78 ± 0.01 0.16 ±0.01 

^ Number of blue / red GC candidates without background corrections and Poisson errors. 

^ Number of blue / red GC candidates background-corrected and Poisson errors including the background errors. 
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Figure 6. Colour distribution functions for the GC candidates around 
NGC 1387 (upper panel), NGC 1379 (middle panel), and NGC 1374 (bot- 
tom panel). Dash-dotted lines show the raw colour distributions and solid 
lines the background-coirected ones. The histograms of the data are plot- 
ted with dashed lines and the background colour distributions with dotted 
lines. Note that the displayed colour distributions and histograms span over 
a colour range slightly larger than the one adopted for the GC selection (i.e. 
0.8 ^ (C-Ti) ^2.3). 

of three. They also reach redder colours as compared to the other 
GCSs. 

From Figs. 5 and 6 and from the data depicted in Table 3, 
the GCS of NGC 1379 shows a similar number of blue GC candi- 
dates as red ones, although suggesting a larger number of red ones 
(40 blue/49 red). The same trend is present for the GC candidates 



around NGC 1374, where the relative number of blue/red clusters 
results 61/73. 

The fractions of blue/red cluster candidates for the three galax- 
ies are 32 per cent for NGC 1 387, 82 per cent for NGC 1 379, and 84 
per cent for NGC 1374. Moreover, all three galaxies present more 
red clusters than blue ones, which is quite peculiar as compared 
to GCSs in general. One may suspect that either some kind of in- 
teraction might have been causing a loss of blue GCs from these 
galaxies, or that an enhanced cluster formation rate was responsi- 
ble of an overabundance of red GCs in any of them, or both. We 
will come back to this point in the discussion. The results of the 
two-Gaussian fits applied to all colour distributions are also given 
in Table 3. The blue peaks agree within the errors and lie in the 
range (C — Ti) = 1.25 — 1.28; it is possible to compare them 
with the blue peak colours determined for GCSs of other Fornax 
galaxies: (C — Ti) = 1.3 2 for GC candidates between r = 1.8 and 
4.5 arcmin in NGC 1399 toirsch et al."2003a*), (B ~ I) = 1.6 for 
GCs with r < 4 arcmin in NGC 1404 (Forbes et al. 1998) which 
corresponds to (C — Ti) = 1.24 a pplying the colour conversions 
proposed bv lForbes & ForteH200lh . and (C - Ti) = 1. 35 for can- 
didate s between r = 0.5 and 5 arcmin in NGC 1427 iForte et all 
I2OOII) . They all agree among each other within 0. 1 mag. 

On the other hand, the red peaks of the three target galax- 
ies show larger differences among each other. Their range is (C — 
Ti) = 1.65 - 1.78, with a bluer peak for NGC 1379. The red peaks 
(within the same radial ranges and obtained by the same authors) 
have (C - Ti) = 1.79 for NGC 1399, (B - /) « 2.1 which cor- 
responds to (C - Ti) = 1.76 for NGC 1404, and (C - Ti) = 
1.76 for NGC 1427. They agree with the NGC 1387 value while 
peaks for NGC 1374 and NGC 1379 are bluer, more noticeable for 
NGC 1379. 



3.2 Spatial distributions 

3.2.1 Aziinuthal distributions 

In order to analyse the radial density distributions of the GCSs, 
it should first be determined if they are spherical or if they show 
some ellipticity. Fig. 7 displays the azimuthal distributions of GCs 
for the three systems with respect to position angle (measured from 
North to East). No evidence for ellipticity can be seen in case of 
the NGC 1387 system. The GCS in NGC 1379 does not show any 
clear azimuthal dependence either within the uncertainties; as this 
galaxy is classified as EO, this is what one expects. 

The same holds true for NGC 1374. It is classified as an El 
galaxy and has a ellipticity smaller than 0.1 as calculated with the 
ELLIPSE task within IRAF. So we expect again sphericity and all 
we can see in Fig. 7 seems to be an excess of GC candidates for 
position angles greater than 180° which is probably due to GCs 
from NGC 1375. 
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Figure 7. Azimuthal distribution of GC candidates in NGC 1387 (upper 
panel), NGC 1379 (middle panel) and NGC 1374 (lower panel). Histograms 
bins span over ^ 17° and Poisson errors are shown for them. 



Th ese results are consistent with those of lKissler-Patig et alJ 
h997j) . 



3.2.2 Radial distributions 

Radial number density profiles of the GCSs are presented in Fig. 8; 
they include the profiles uncorrected for background contamination 
as well as the background-corrected profiles. For the latter case, 
blue and red candidates are shown separately and are tabulated in 
Table 4. In all cases, the Poisson errors include the uncertainties of 
the raw counts and of the background. The T\ luminosity profiles 
of the corresponding galaxies, in a proper scaling, are also shown as 
well as the densities of background objects determined in Section 
2.2. We remind that the adopted cut-off magnitude is T\ = 23.7. 

The innermost bins, which correspond to radial distances 100 
< r < 200 pix (0.45 - 0.68 arcmin) from the galaxy centres, are 
affected by the light of the galaxies so we apply an 80 per cent 
completeness factor to correct for this 'geometric' incompleteness. 

The density profile of NGC 1387 clusters shows that the raw 
counts profile flattens at a level that agrees with the background 
density estimated above. The background-corrected density profile 
was fit by a power law of the form A'^ = a.r^ in the range 54 arcsec 
to 162 arcsec (2nd. to 5th. bins in Table 4), excluding the innermost 
bin and the outer regions where there are very few candidates left. 
The slope obtained for all GC candidates was b — — 2.4±0.3, while 
for blue and red candidates separately they were b = —2.0 ± 0.5 
(blue clusters) and b=~2.6±0.3 (red clusters). The GC density 
profile in NGC 1387 is clearly dominated by the red clusters, which 
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Figure 8. Radial density profiles for GC candidates in the field of 
NGC 1387 (upper panel), NGC 1379 (middle panel), and NGC 1374 (lower 
panel). Open circles show the radial distributions uncorrected for back- 
ground contamination and filled circles the background-corrected distri- 
butions, which are split in blue candidates (open squares) and red candi- 
dates (open triangles). Solid lines represent the Ti brightness profiles of the 
galaxies multiplied by —0.4 and shifted by an arbitrary constant, dashed 
lines are the power-law fits to the GC radial density profiles, dot-dashed 
lines are the power-law fits to the galaxy brightness profiles, and dotted 
lines indicate the background densities. The innermost bins have been cor- 
rected by geometric incompleteness. 



show a steeper profile than the blue candidates, a result in com- 
mon with other GCSs as, for example, the one s in the already men - 
tioned Forn ax early type galax ies NGC 1399 jDirecl^^lj2003£ ), 
NGC 1404 fror besetalll9 9^. and NGC 1427 l lForte et all200ll) . 
The three profiles (for all, red and blue candidates) are also fit rea- 
sonably well by a law logN = c + d.r^^^ , within the same range 
used for the power law. 

The density profiles of the GC population in NGC 1379 look 
similar to the ones for NGC 1387, except that no dominant subpop- 
ulation is present in this case as red candidates are only slightly 
more abundant than blue ones. The raw counts profile levels out 
at the density value corresponding to the background objects, es- 
timated near this galaxy. The background-corrected profile can be 
fit by a power law of slope b — — 2.6±0.5in the same radial 
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Table 4. Radial number densities of blue and red GC candidates (in arcmin ^, and background-corrected) in the galaxies NGC1374, NGC 1379 and 
NGC 1387. 



r [arcmin] 


Tvr blue 
^^NGC 1374 


M red 

^^NGC 1374 


]\T blue 
^^NGC 1379 


^^NGC 1379 


]vr blue 
^^NGC 1387 


^^NGC 1387 


0.45 


-0.90 


9.92 ±2.78 


16.34 ±3.34 


12.69 ±3.00 


11.98 ±2.86 


6.49 ± 2.36 


18.22 ±3.53 


0.90 


- 1.35 


4.41 ± 1.41 


8.73 ± 1.72 


3.02 ± 1.14 


5.19 ± 1.33 


2.06 ± 1.14 


8.66 ± 1.72 


1.35 


- 1.80 


2.61 ± 1.01 


2.22 ±0.81 


0.73 ±0.64 


1.33 ±0.64 


1.57 ±0.90 


5.29 ± 1.17 


1.80 


-2.25 


1.09 ±0.72 


1.40 ±0.61 


1.03 ±0.61 


0.76 ± 0.46 


0.94 ±0.72 


2.03 ± 0.70 


2.25 


-2.70 


0.0 ±0.50 


0.00 ±0.25 


0.36 ± 0.46 


0.54 ±0.38 


0.40 ±0.59 


1.23 ±0.54 


2.70 


-3.60 


0.20 ±0.35 


0.14 ±0.22 


0.45 ± 0.30 


0.15 ±0.19 


0.0 ±0.35 


0.13 ±0.23 


3.60 


-4.50 


0.0 ±0.26 


0.30 ±0.22 


0.07 ±0.23 


0.0 ±0.13 


0.0 ±0.28 


0.0 ±0.19 


4.50 


-5.85 


0.0 ±0.24 


0.12 ±0.16 


0.01 ±0.17 


0.01 ±0.11 


0.21 ±0.24 


0.0 ±0.14 


5.85 


- 7.20 


0.0 ±0.18 


0.0 ±0.14 


0.16 ±0.16 


0.08 ±0.11 


0.45 ±0.22 


0.03 ± 0.13 



range as in NGC 1387 (2nd. to 5th. bins in Table 4); from the fits 
of blue and red candidates we obtain b = —2.3 ± 0.8 (blues) and 
h = — 2.9±0.4 (reds). The red clusters seem to be again more con- 
centrated than the blue clusters, although the uncertainty is rather 
large. As in the case of NGC 1387, an r^^'^ law provides a good 
fit for the density profiles of all candidates together as well as for 
blues and reds separately. 

Finally, the GC density profiles for NGC 1374 present a dif- 
ferent behaviour, as depicted in Fig. 8. The raw counts profile does 
not flatten at the density value of the nearby estimated background 
objects due to the presence of GCs which probably belong to neigh- 
bouring galaxies; the candidates from NGC 1375 may be the re- 
sponsible of the excess at r > 162 arcsec or 600 pix (see Table 4) 
while some candidates from NGC 1373 may be also present at the 
two outer bins. To avoid this extra contamination we restrict the 
fit of the background-corrected counts by a power law to a smaller 
radial range 54 arcsec < r < 135 arcsec (2nd. to 4th. bins in Ta- 
ble 4), from which we obtain slopes h — —2.3 ± 0.2 (all candi- 
dates), 6 = -2.3 ± 0.6 (blue clusters) and b = -3.2 ± 0.6 (red 
clusters). 

Th e slope co e fficient s for the power law fits found by 
fKissler-P atig et alJ <1997j) . for all candidates, in NGC 1387, 
NGC 1379 and NGC 1374 are: -2.2, -2.1 and -1.8, respectively; 
they agree within the uncertainties with our values except for 
NGC 1374 for which they obtain a shallower slope, probably due 
to the influence of the nearby contaminating GCSs. 

It is difficult to compare our results for NGC 1379 with the 
HST based ones from lElson et alJ ^99^ due to the different spa- 
tial coverage, limiting magnitude, and comparison field (see Sec- 
tion 3.1). Elson et al. argue that the GCS is lost in the back- 
ground at 80 arcsec while, on the other side, they are able to anal- 
yse its behaviour even within 10 arcsec; instead, our density pro- 
file suggests that the GC density reaches the background level at 
about 200 arcsec whi le we cannot study the GCs inside 27 arcsec. 
ICapuzzo-Dolcetta"&^ Donnarumma ( 2001) assign the difference be- 
tween the surface brightness profile of NGC 1379 and its GCS pro- 
file in the inner region (r < 50 arcsec) to a loss of GCs caused 
by dynamical friction and tidal interaction with the galactic nu- 
cleus; in our case the innermost bin in the radial density profile of 
NGC 1379, corrected for geometric incompleteness, does not show 
any noticeable difference with the galaxy light profile. 

For the three galaxies under study, the radial density pro- 
files GC candidates follow the Ti light profiles of the respective 
galaxies (at least out to the radial distance at which they become 
undistinguishable f r om the background), as already pointed out by 
lKissler--Patig et alj Jl997iJ) . The power law indices obtained for 



them are b = -2.2 ± 0.1 for NGC 1387, b = -2.7 ± 0.1 for 
NGC 1379, and b = -2.4 ± 0.1 for NGC 1374. 



3.3 Luminosity functions 

The luminosity functions for the GC candidates (GCLFs) in 
NGC 1387, NGC 1379 and NGC 1374 are presented in Fig. 9. As 
the images used in this study are not very deep, and the complete- 
ness sets the faint magnitude limit at Ti — 23.7, we are not able 
to reach the turn-over magnitudes (TOMs) of the GCLFs with the 
required precision to derive distances. Instead, we opt to investigate 
whether these TOMs show the universal behaviour. 

The GCLFs are based on the magnitude distributions corrected 
for photometric incompleteness and contamination by the back- 
ground with bins of 0.6 mag for NGC 1387 and NGC 1379 and 
bins of 0.7 mag for NGC 1374 for which we choose a slightly larger 
bin to diminish the errors as incompleteness seems to be affecting 
the fainter bins more in NGC 1374 than in the other galaxies (see 
Fig. 9). As said above, an inner limit of r = 27 arcsec is considered 
for all three GCSs and the geometric incompleteness in the 100 - 
200 pix bins have not been corrected. 

As we are unable to construct the GCLFs beyond the TOMs, 
it is not possible to fit Gaussians leaving both mean value and dis- 
persion as free parameters. We rather prefer to fix the dispersion 
at CT — 1.2, a val ue already obtained on t he basis of Washing- 
ton photometry by 'Ostro v et all h998 f> and 'Pi rsch et all i2003j) 
for NGC 1399, and by For teetalJ J2001) for NGC 1427, and just 
estimate indicative TOMs by means of the Gaussian fits. 

The fit to all the NGC 1387 GC candidates gives an indica- 
tive TOM at Ti = 23.54 ± 0.06, and for the red subpopula- 
tion at Ti = 23.53 ± 0.18, in agreement with the former value 
(more than 75 per cent of the GCs belong to the red subpopu- 
lation in this galaxy). We can also estimate the ex pected TOM 
if we adopt the SBF distance modulus obtained by Ijensen et alJ 
(200i) for this galaxy (m - M) = 31.38 ± 0.26 and ac- 
cept a universal TOM. The absolute TOM estimated bv * Richtie3 
(2003) as a weighted avera ge of the TOMs for the Milky Way and 
M31 llHarrij20C)li;iBarmbv. Huchra & BrodiH200lL respectively) 
is Mvq = —7.46 ± 0.18. By means of the mean colour of the 
whole GC sample (C — Ti ) = 1.5 4 ± 0.34 and the colour conver- 
sions from iForbes & Fort3 i200ll) . we can transform the V band 
TOM into an R band one Mr = —8.00. Thus, the expecte d appar- 
ent TO M is _R = 23.38 (which corresponds to Ti = 23.40 JCeisleil 
Il99d) ); our observed value is 0.14 mag fainter. 

In the case of NGC 1379 and NGC 1374, we perform fits to 
the whole population of GC candidates only, without separating 
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Table 5. Total GC populations Nqc, integrated absolute magnitudes My 
and specific frequencies Sn for the three target galaxies. 



Figure 9. Ti (dereddened) luminosity functions for the GC candidates 
in the field of NGC 1387 (upper panel), NGC 1379 (middle panel), and 
NGC 1374 (lower panel), for r > lOOpix (27 arcsec), and corrected for 
contamination by the background and photometric incompleteness. SoHd 
lines show Gaussian fits to all candidates, and the dotted line the Gaussian 
fit only to the red candidates (open triangles); all fits were peii'ormed with 
a fixed dispersion ct = 1.2 and for magnitudes up to Ti = 23.7, according 
to the 50 per cent completeness limit (in the shaded region the results are 
uncertain). Geometric incompleteness in the 100 - 200 pix annuli have not 
been corrected. 



between blue GCs and red GCs. There are no dominant subpopu- 
lations and the number of only blue or red candidates would be too 
low. The indicative TOMs then result in Ti = 23.33 ± 0.14 for 
NGC 1379, and Ti = 23 07 ± 0.11 for NGC 1374. The SBF dis- 
tances from lJensen et alj<2003h are (m - M) = 31.35 ± 0.15 for 
NGC 1379 and (m - M) = 31.32 ± 0.13 for NGC 1374. Follow- 
ing the same procedure, and by means of the mean colour of the GC 
populations: (C-Ti) = 1.42±0.34 forNGC 1379 and (C-Ti) = 
1.50 ± 0.28 for NGC 1374, we estimate the expected absolute 
TOMs as Mr = -7.97 and Mr = -7.99, and the corresponding 
apparent TOMs as Ti = 23.40 and Ti = 23.35, forNGC 1379 and 
NGC 1374, respectively. For NGC 1379 the observed and expected 
TOMs agree within the errors while for NGC 1374 the observed 
TOM is about 0.3 mag brigh ter. The observed TOMs for NGC 1374 
given bv lKohle et ai]|l996l) are V = 23.52 from a Gaussian fit and 



Galaxy 


Ngc 


Mv 




Sn 


NGC 1374 


360 ± 17 


-20.4 ± 


0.1 


2.4 ±0.5 


NGC 1379 


225 ± 23 


-20.6 ± 


0.2 


1.4 ±0.4 


NGC 1387 


390 ± 27 


-20.9 ± 


0.3 


1.8 ±0.7 



V = 23.44 from a ts-function fit. The expected TOM according to 
the SBF distance isV = 23.84, so in this case again the observed 
TOMs turn out to be 0.3 and 0.4 mag brighter, respectively, than 
the expected one. 



3.4 Total populations and specific frequencies 

In order to calculate the specific frequencies we need the total GC 
populations and the absolute V magnitudes of the galaxies. The to- 
tal GC populations can be estimated from the luminosity functions, 
if we accept that they are fully Gaussian and the corresponding 
mean values and adopted dispersions are known. However, for all 
the GCSs under study, we lack information on the central 100 pix 
area, and also the surface density values for the first radial bin (100 
pix - 200 pix) are affected by geometric incompleteness that was 
not corrected in the GCLFs. We will then estimate, for each sys- 
tem, the population within the inner 100 pix assuming a constant 
surface density determined from the power law fit at r = 100 pix; 
we Justify this assumption for the three GCSs following lElson et alJ 
(1998) who found for NGC 1379 that the surface density profile of 
the GCs flattens out inwards r = 30 arcsec, which corresponds to 
111 pix on our scale. In addition to this, a correction that should 
be taken into account is the difference between the observed radial 
density obtained for the first bin (uncorrected for the incomplete- 
ness) and the calculated radial density from the power law fit for 
the same bin. Thus, the inner 100 pix population and the correction 
for the 100 - 200 pix range will be added, for each GCS, to the 
populations obtained from the GCLFs. 

First, we calculate the raw number of GCs in the bright halves 
of the Gaussians (taking the TOM as upper magnitude limit), and 
correct them for contamination by the background and photometric 
incompleteness; these counts are then doubled to cover the whole 
Gaussian. The errors are calculated on the basis of the Poisson 
uncertainties of the raw and background counts, and the effect 
of the incompleteness. In this way, we obtain 274 ± 18 GCs for 
NGC 1387, 160 ± 11 GCs for NGC 1379, and 237 ± 12 GCs for 
NGC 1374. Secondly, the fraction of the GC populations contained 
in the inner 100 pix corresponds to 72 ± 15 GCs for NGC 1387, 
50 ± 15 GCs for NGC 1379, and 76 ± 8 GCs for NGC 1374, where 
the errors take into account the errors of the fit. Finally, the cor- 
rections to be added due to the differences between observed and 
calculated densities in the first radial bins of each system, and their 
corresponding errors, are of 44 ± 14 GCs for NGC 1387, 15 ± 13 
GCs for NGC 1379, and 47 ± 9 GCs for NGC 1374. Thus, the to- 
tal estimated populations result in 390 ± 27 GCs for NGC 1387, 
225 ± 23 for NGC 1379, and 360 ± 17 GCs for NGC 1374, where 
the uncertainties are estimated as a combination of the uncertainties 
in the thre e steps just de scribed; they are shown in Table 5. 

Kissl er-Patig et alj {l997a) obtained 389 ±110 GCs for 
NGC 1387 314 ± 63 GCs for NGC 1379, 410 ± 82 GCs for 
NGC 1374. lElson et"ai]jl998h quoted 436±30 GCs for NGC 1379. 
Our results for NGC 1387 and NGC 1374 are in good agreement 
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Figure 10. R integrated magnitudes vs. radius for NGC 1374. The solid 
line shows our results and solid squares plot the aperture photometry from 
[Poulain & Nieto, U994 ). The vertical dashed line indicates the limiting ra- 
dius where R integrated magnitude has been measured for further calcula- 
tions. 

with lKissler-Patig et alJh997ah . particularly for NGC 1387. In the 
case of NGC 1379, our estimated GC population agrees barely 
within the errors with the one from Kissler-Patig et al., and it is 
about 50 per cent below the estimation from Elson et al. This large 
difference is likely a consequence of the different background: as 
the Elson et al. background is far away from NGC 1399, it lacks the 
contribution of its CGS and this might turn into a overabundant GC 
population. 

The galaxies Ti brightness profiles were integrated out to the 
respective radii adopted for each GCS and yield T\ = 9.98 ± 0.06 
for NGC 1387 (r < 3.1 arcmin), Ti = 10.29±0.05 for NGC 1379 
(r < 2.7 arcmin), and Ti = 10.39 ± 0.04 for NGC 1374 (r < 2.2 
arcmin). The uncertainties of the integrated Ti magnitudes arise 
basically from the standard errors of the respective sky determi- 
nations. Our Ti images are saturated at the three galaxy centres, 
over small regions with radii r ~ 40 pix (about 0.2 arcmin). As 
we measure integrated magnitudes at larger limiting radii (r = 2.2, 
2.7 and 3. 1 arcmin), the proportion of light missed at the centres is 
negligible in comparison to the light integrated up to such limits. 
The Ti integrated magnitudes were transformed to R magnitudes 
through the zero-point difference (R — Ti ~ —0.02) with the pur- 
pose of testing them against published R multi-aperture photome- 
try from Poulain (1988.) and .Poulain & Nieto 1.1994.) . In Fig. 10 we 
plot published and calculated integrated magnitudes for NGC 1374, 
the galaxy with the smallest limiting radius, to show that the dif- 
ferences are evident at inner radii but tend to vanish for r > 0.75 
arcmin, that is, at a radius much shorter than the limiting one. We 
do not include the same plots for the other two galaxies as they look 
very similar. 

Adopting as representative colours for the galactic haloes 
of NGC 1387, NGC 1379 and NGC 1374 the values obtained 
from the surface brightness profiles: (C — Ti) = 1.61, 1.49, 
and 1.45, respectively , and applying the colour conversions from 
iForbes & Fortd i200lh . we calculate the corresponding V magni- 
tudes integrated up to those limiting radii asV = 10.52, 10.80, and 
10. 89, respectively. B y means of the adopted SBF distance mod- 
uli i Jensen et all2003l) . we finally obtain the integrated absolute V 
magnitude (up to the limiting radii) for NGC 1387, NGC 1379 and 
NGC 1374: Mv = -20.86± 0.3, -20.55± 0.2, and -20.43± 
0.1, respectively, whose uncertainties include the uncertainties of 



the apparent integrated magnitudes and the distance moduli. In this 
way, the specific frequencies calculated with the GC populations 
and absolute V magnitudes, both estimated up to the same radii, are 
Sjv = 1.8 ± 0.7 for NGC 1387, Sn = 1.4 ± 0.4 for NGC 1379, 
and Sn ~ 2.4 ± 0.5 for NGC 1374; the uncertainties are calcu- 
lated on the basis of the uncertainties of the absolute magnitudes 
and of the GC populations. Absolute V magnitudes and specific 
frequencies are also given in Table 5. 



4 DISCUSSION 

4.1 Testing the tidal stripping scenario 

The three galaxies under study present several characteristics that 
point to a particular interpretation, (i) Harris (2003) quote a 'typi- 
cal' specific frequency of Sjv ~ 4 for elliptical galaxies but the val- 
ues determined for NGC 1379, and NGC 1374 are distinctly smaller 
(NGC 1387 is an SO galaxy, where a small S'jv -value is not un- 
usual) (ii) The GCSs of these three galaxies show a rather small 
proportion of blue GCs with respect to the total GC population: 
24 per cent for NGC 1387, 45 per cent for NGC 1379, and 43 per 
cent for NGC 1374. (iii) The radial distributions of blue clusters are 
less concentrated than the red ones with respect to the galaxy light 
which implies a higher probability for the blues to get lost during 
tidal stripping processes, (iv) The galaxies are all located close to 
NGC 1399, a giant elliptical galaxy with a very rich GCS. These 
properties make it likely that these three galaxies have experienced 
tidal stripping of their blue GCs by NGC 1399. Kissler-Patig et ^ 
( 1999) already suggested that the rich cluster system of NGC 1399 
have partly been formed by tidal stripping of GCs from neighbour- 
ing g alaxies. 

Bek ki etaP i200 3') perform numerical simulations on the 
tidal stripping and accretion of GCs for the case NGC 1404 and 
NGC 1399 in Fornax. They propose to test the influence of tidal 
stripping in low Sn ellipticals that might be involved in such pro- 
cess, looking for 'tidal streams' of intracluster GCs along the orbit 
of their former parent galaxy, and checking if there is some corre- 
lation between the distance of an elliptical galaxy from the center 
of the cluster and Sn- With regard to the first test, the GCs around 
NGC 1399 present asymmetries in their azimuthal projected distri- 
bution and a 'bridge' of blue GCs in the direction to NGC 1387 
appear quite distinctly, pointing to an overdensity of GCs in the 
particular direction. We cannot assure that these GCs have been in 
fact lost by NGC 1387 as its orbit around NGC 1399 is not known, 
but being the blue globulars less bound to NGC 1387 than the reds, 
due to their larger mean distance from the host galaxy, they have a 
higher probability of being lost. 

The other test proposed by iBekki et'aP j2003) is to search 
for a correlation between the Sn of the elliptical galaxies and 
their projected distances to NGC 1399. Fig. 11 shows such rela- 
tion for the NGC 1374 and NGC 1379 as well as for NGC 1404 
and NGC 1427. Other Fornax galaxies close to NGC 1399 like 
NGC 1380 and NGC 1387 are not included as they are lenticulars. 
There seems to be a trend in the sense that galaxies with larger pro- 
jected distances from the cluster center tend to have higher specific 
frequencies, but there are only four points in t he graph. Previous re - 
sults related to this test have been presented bv lForbes et alJi^997^ . 

Further evidence for this may be expected from kinematical 
studies of the NGC 1399 system at radii larger than presently avail- 
able velocities <Richtler et all2004l) . 
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Figure 11. Specific frequency vs. projected distance from NGC 1399. From 
left to right, the data are from NGC1404, NGC 1379, NGC 1374 and 
NGC 1427. The S n were taken from jorbes et aU tl998i) . this paper and 
lForteetalJ<200lh . respectively. 



Table 6. Total visual magnitudes and (C — Ti ) colours of the blue and red 
GC peaks for galaxies with Washington photometry. 



Galaxy m — M' 



(C-Ti) 



'red pk 



Sour.'^ 



NGC 1374 


31.32 


-20.22 


1.27 ±0.02 


1.70 ±0.03 


1 


NGC 1379 


31.35 


-20.36 


1.28 ±0.02 


1.65 ± 0.01 


1 


NGC 1387 


31.38 


-20.66 


1.25 ±0.01 


1.78 ± 0.01 


1 


NGC 1399 


31.02 


-21.53 


1.32 ±0.05 


1.79 ± 0.03 


2 


NGC 1427 


31.02 


-20.11 


1.35 ±0.07 


1.76 ± 0.07 


3 


NGC 3258 


32.53 


-21.23 


1.36 ±0.04 


1.70 ±0.02 


4 


NGC 3268 


32.71 


-21.41 


1.31 ±0.02 


1.66 ±0.03 


4 


NGC 3923 


31.80 


-22.11 


1.47 ±0.025 


1.87 ±0.025 


5 


NGC 4472 


31.20 


-22.82 


1.31 ±0.04 


1.81 ±0.04 


6 


NGC 4636 


31.24 


-21.73 


1.28 ±0.02 


1.77 ± 0.02 


7 



^ Distance moduli from GC luminosity functions (when available) or SBF 
distances i jensen et a l. 1999; Tonrv et al. 2001; Jensen et al. 2003). 
^ My calculated with V total magnitudes from RC3 catalog and the dis- 
tance moduli quoted in the second column. 

— Ti) peak co lour s sources. 1: t his pap er, 2: [Pirsch et alJ j 2Q03ah . 
3: Forte etal. (2001), 4: Dirsch et al. 1 2003bl). 5: IZepfetalJ Il995t) . 6: 



"leisleretaU 11996.1 . 7: .Dirsch et aL^200.5il 



4.2 The peculiar colour distribution of the GCS in NGC 1387 

The colour distribution of the NGC 1387 GCS (see Fig. 6) is pecu- 
liar due to the large number of red clusters and also because they 
extend to redder colours than in most GCSs. One possibility to ex- 
plain such distribution is the presence of background galaxies that 
might have been erroneously classified as red GC candidates. No 
contamination of this type is expected from the results of the back- 
ground adopted (region 'a' in Fig. 4); it can be seen in Fig. 6 that 
the contribution of objects redder than (C — Ti) — 1.45 is very 
low and does not affect the shape of the red GCs colour distribu- 
tion. Even if such unresolved galaxies were present in the same 
area as the GCS, it is shown in Paper I that they would have con- 
taminated the GCs colour range in — / but not in (C — Ti) as 
they appear at (C — Ti) < 0.8 (see Paper I, their fig. 7). Objects 
with (C — Ti) < 0.8 and Ti > 23 do not show any concentra- 
tion towards the target galaxies, like the GC candidates, but they 
present an homogeneous projected distribution over the whole field 
considered in the present paper. Thus, we expect the contamination 
of the red candidates by background galaxies to be low. 

The lack of blue GCs may be explained by tidal stripping pro- 
cesses, as mentioned above, but NGC 1387 is an SO galaxy and 
alternative explanations may be valid as well. The GCS of the 
SO galaxy NGC 1380 in the Fornax cluster has been studied by 
[Kissler-Patia et al. ( 1997b) and it presents several properties sim- 
ilar to the ones we describe in the NGC 1387 system. The colour 
distribution of NGC 1380 GCs is clearly bimodal and their {B~V) 
and {B — R) colour histograms (their figs. 3 and 4) show a large 
number of red GCs that extend to redder colours than the metal-rich 
population in the Milky Way, like in NGC 1387. The fraction of 
blue clusters over the total population is small (32 per cent) as well 
as the specific frequency obtained by Kissler-Patig et al. Sn = 
1.5 ± 0.5. These similarities between the GCSs of NGC 1380 and 
NGC 1387 are also evidence that the peculiar GCS colour distribu- 
tion in NGC 1387 is hardly due to background galaxies. 

On the other hand, the blue population of NGC 1380 shows 
a spherical distribution around the galaxy while the red one fol- 
lows the stellar light in ellipticity and position angle so they can 
be associated with the galactic bulge and disc. When we study sep- 



arately the radial distributions of blue and red GCs in NGC 1387, 
the blues show no sign of ellipticity and the reds only a marginal 
evidence, too weak to be significant. In this way, we cannot dis- 
tinguish between a bulge and disk population. The formation of 
SO galaxies is still a matter of discussion. As a class, they ar e 
probably not simply gas-stripped spirals (Bu rstein et alJ EoO^ . 
However, interaction between galaxies is likely to play a role in 
their formation history (e.g. Bicker et al. 2003; Shiov a et all2004l 
iFalcon-Barroso. Peletier & Emsellem et aL. 2004). This interaction 
could have triggered a starburst in NGC 1387, causing the abun- 
dant formation of metal-rich clusters. Since NGC 1387 is seen al- 
most face-on, a cluster disk population is not distinguishable from a 
halo population by its projected spatial distribution. A radial veloc- 
ity study could reveal whether the metal-rich clusters in NGC 1387 
show a disk kinematical behaviour. 



4.3 The GCS colour distributions compared to the 
luminosities of the parent galaxies 

In order to analyse the trend of the (C — Ti) colours of the blue 
and red peaks of GCSs with luminosity of the corresponding par- 
ent galaxies we present in Table 6 the data obtained from previous 
Washington photometric studies and from this paper. The sources 
are detailed in the table. 

Blue and red peak colours versus visual absolute magnitudes 
are plotted in Fig. 12. By means of linear fits performed to the 
data we obtain a slope of — 0.02± 0.02 for the blue peaks and 
— 0.04± 0.02 for the red ones. The general appearance resem- 
bles that found in {V — I) i Larsen et al. 2001 ; Kundu & Whitmo^ 
12001 Jd) : the colour of the red peak depends slightly on the lumi- 
nosity of the host galaxy, while the blue peak shows a much shal- 
lower dependence, reasonably consistent with being constant. The 
scatter is of the order which can be expected for a inhomogeneous 
sample and small number statistics cannot be ignored either. The 
only previous Washington photometric study that has not been in- 
cluded in the sample of GCSs depicted in Table 6 is NGC 33 11 
as it seems to pose a problem. ISecker et alj il995h claimed that 
NGC 33 11 lacks the population of metal-poor clusters. However, 
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Figure 12. Blue and red (C — Ti) peak colours of GC colour distribu- 
tions in galaxies with Washington photometiy vs. luminosity of the parent 
galaxies. The solid lines show the linear fits to the red and blue GC peaks, 
respectively. 

iBrodie. Larsen & Kissler-Padd \2QQ& i found the (V - 7) colour 
distribution quite normal. Since the photometric quality of the 
Washington photometry was doubtful (Geisler, priv. communica- 
tion) we excluded this GCS of the selected sample. 

That the red peaks have a slight dependence on the host 
galaxy's luminosity seems to be plausible: brighter galaxies possess 
more metal-rich clusters. The well-known relation between bright- 
ness and metallicity for early-type galaxies should be somehow re- 
flected in the metallicities of the metal-rich clusters. T he constancy 
of the blue peak has been a matter of discussion. iLarsen et alj 
^20011) find that the blue peak becomes bluer with decreasing 
galaxy luminosity (by an amount of 0.016 mag in {V — I) per B- 
magnitude of the host galaxy), while Kundu & Whitmore ( 2001a b) 
do not see this effect. Given the superior metallicity resolution of 
the Washington system over [V — I), we would expect to uncover 
a relation if there was any. On the other hand, a constant blue peak 
does not necessarily mean that all metal-poor GCs have the same 
metallicity distribution. As mentioned in Dirsch et al. (2003a), a 
non-linear colour-metallicity relation in combination with a certain 
scatter around the mean relation (by second parameters or photo- 
metric errors) can produce a quite universal blue peak, if a given 
metallicity distribution is rebinned in colour. The c olour-metallicity 
relati on for the Washington system is non-linear iHarris & Harris! 
l2002h . Due to its large scatter, the {V — I) colour-metallicity re- 
lation is difficult to calibrate in the metal-poor regime of Galactic 
globular clusters, but it is reasonable to assume that it also becomes 
non-linear. 

This does not n e cessar ily contradict the results of 
IStrader. Brodie & Forbes! i2004) . who found a significant re- 
lation between the mean blue colours in {V — I) and galaxy 
luminosity. They include very faint dwarf galaxies whose cluster 
systems contain only a few clusters and where the quotation of a 
'peak' is meaningless. 



5 CONCLUSIONS 

The colour distributions of the GCSs in the low-luminosity galax- 
ies NGC 1387, NGC 1379, and NGC 1374 show a distinctly bi- 



modal appearance; these bimodalities have not been detected in 
previous work. NGC 1387 even shows a colour bimodality more 
notably than known from any other GCS. That bimodality is not 
an exclusive property of the GCSs of giant elliptical galaxies, 
has already been shown in an other Fornax low-luminosity galaxy: 
NGC 1427 <Forte et aljboOl h. Blue and red GC candidates also 
show different spatial distributions, the red GCs being more cen- 
trally concentrated than the blue ones. The appearance of bimodal- 
ity in our Washington colours probably is a consequence of the 
superior metallicity resolution of the Washington system with re- 
spect to other broad band colours. In fact, all GCSs observed so 
far in the Washington system revealed bimodal distributions, while 
in {V — I), perhaps due to the inferior metallicity resolution, a 
certain fraction remained unimodal or unclear (e.g. !Larsen et alJ 
1200 it iKundu & Whitmore! EoOlalbh . One may suspect that a bi- 
modal colour distribution is a property of all GCSs, as long as there 
is a statistically significant to use this concept. It is therefore desir- 
able to investigate a larger sample of low-luminosity galaxies in the 
Washington system. 

With the aid of the luminosity functions we estimate, for these 
three low-luminosity galaxies, GC populations between 200 and 
400 GCs and specific frequencies in the range Sn = 1.4 - 2.4 , 
i.e. sm aller than the typical values for elliptical galaxies I Harris 
!2003t) . though we remind here that NGC 1387 is a SO galaxy, not 
an elliptical. Though it is likely that the GCSs under study owe 
the characteristics of their GC colour distributions and low specific 
frequencies to some interaction processes, as discussed in the pre- 
vious section, it is worth to analyze them, in the following, under 
the hypothesis that they are just regular systems whose properties 
have not been affected by, for instance, tidal stripping. 

The specific frequencies provide us with valuable information 
related to the efficiency of the galaxies at forming GCs as well as 
the role the environment might have played during the evolution. 
For instance. jForbes et alJi200ill found that specific frequencies of 
red GCs estimated with respect to (only) the 'bulge' luminosity of 
spiral galaxies Sjv = 0.5 - 1.5 are similar to the ones of red GCs in 
field elli pticals estimat ed in the usual way Sn = 1-3 (see refer- 
ences in !Forbes et alj2 00 1 ) . It is interesting to note that the specific 
frequencies we obtain for the red candidates around the ellipticals 
NGC 1374 (Sn"^ = 1.4 ± 0.3) and NGC 1379 (S^"^ = 0.7 ± 0.3) 
are also close to these Sn ranges which can be considered as an 
evidence in favour of the Forbes et al. claim that metal-rich GCs in 
spirals and ellipticals may have the same origin: they are formed 
with the bulge stars, assuming that ellipticals are bulge-do minated 
galaxies. The numerical simulations performed by iBeaslev et alJ 
( 2002) on the formation of GCs within a hierarchical Universe also 
provide theoretical support to this statement. 

Let us turn no w to the predictions particularly related to low- 
luminosity galaxies. lForbes et alJil997h discussed different models 
for the origin of GCs in giant elliptical galaxies and suggested, as 
the most probable scenario, the formation of GCs 'in situ' in two 
star formation episodes. One of the predictions that results from 
the Forbes et al. analysis and can be tested through observations, is 
that low-luminosity, low Sn galaxies should have more metal-rich 
than metal-poor GCs. This trend is not follow ed, for example, by 
NGC 1427 (Mv = -20.05. !Forte et all200ll) with a blue/red ratio 
of GCs in the range 3.9 - 4.4, or by NGC 3379 in the Leo I Group 
(Mv = "20.9, Rhode & Zeol 2004) with a 70 % of blue GCs 
with respect to the total globular population. Instead, the three low- 
luminosity galaxies under study present blue/red ratio in the range 
0.3 - 0.8, in agreement with the Forbes at al. predictions. Anyway, 
we should keep in mind that the low number of blues in these GCSs 
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may be originated in another processes, like tidal stripping, and that 
one of the galaxies is of SO type. Thus, our blue/red GC ratios do 
not provide much evidence in favour of the proposed model. 

According to our results, the 'classic' interpretation of bi- 
modality as a consequence of merger events which would add a 
metal-rich GC population to a pre-existing metal-poor GC popula- 
tions, is doubtful in the case of low-luminosity ellipticals. How- 
ever, the hierarchical scenario proposed through simulations by 
iBeaslev et al. (2002), with blue GCs formed in protogalactic cold- 
gas fragments and red GCs formed in subsequent gas-rich mergers, 
leads to bimodal colour distributions in galaxies of all luminosities, 
even low-luminosity ones. According to Beasley et al. predictions, 
the metallicity distributions of blue GCs are very similar in low 
and in more luminous ellipticals. The red GCs in low-luminosity 
ellipticals, as these galaxies suffer fewer mergers followed by star 
formation than the luminous ellipticals, present a larger scatter in 
their mean colours. Both predictions are verified here by the GCSs 
of the ellipticals NGC 1374 and NGC 1379. 

After checking the properties of the GCSs around NGC 1374, 
NGC 1379 and NGC 1387 against previous observational and the- 
oretical results, we suggest that a plausible scenario of formation 
and evolution of GCSs which shows reasonably agreement with 
our result s would be a combin ation of formation in a hierarchical 
Universe jBeaslev et id]|200gl^ and further tida l stripping, or other 
sort of interaction processes " Iporbes et aljl99l for instance), orig- 
inated in their proximity to the giant elliptical NGC 1399. We reit- 
erate that kinematical studies of these GCSs would help to obtain 
more evidence on the proposed statements. 
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